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Abstract 

We discuss the possibility of the existence of a long-lived top 
squark (r <C 10 10 years) and its motivation. If the stop is indeed 
metastable, it forms hadrons. We study properties of the low-energy 
stop-containing hadrons and their signatures in collider experiments. 
PACS: 12.60.J, 14.80.Ly, 14.20 



1 Introduction 



A feature of most discussions on the possibility of observing supersym- 
metry in hadronic physics is that attention has been paid almost exclusively 
to "fragmentation physics." That is predictions have been made of what 
distinctive signature might occur in hadronic jets due to the production and 
subsequent decay of sparticles. This appears due to the expectation that 
the mass of the lightest squarks will be substantially greater than that of 
the lightest super-partners of standard model particles that do not carry 
color. This is the pattern within the conventional standard model and is 
so expected for its presumptive supersymmetric successor. Under this cir- 
cumstance, bound states involving squarks and ordinary quarks would not 
undergo "hadronization," superpartner hadrons ("shadrons") would never 
form and hadron spectroscopy would be irrelevant. 

Although unconventional, we can ask a simple question, "If the mass of 
squarks were comparable to that of the LSP (lightest supersymmetric par- 
ticle), what type of hadrons might be expected?" Part of what prompts 
us to ask this unconventional question has been our formal investigations 
of possible structures that give naive supersymmetric extensions of the 
effective field theory describing low-energy pion physics. These results seem 
to indicate that supersymmetric extensions of these are formally possible. 

The possibility of the existence of supersymmetric hadrons was first 
pointed out by Farrar and Fayet @. Historically, models containing light 
gluinos rather than squarks appeared to be better motivated. Subsequently, 
the gluino containing hadrons ("R- hadrons") have been studied to a great 
extent 0. On the other hand, another option - the squark containing hadrons 
- has not drawn similar attention. We believe that, in the ongoing search 
for supersymmetry, one cannot afford the luxury of ignoring this possibility 
however unlikely it may seem. In this letter, we attempt, at least partially, 
to fill this gap in the phenomenology of supersymmetry. 

2 Metastable Stop 

In this paper we study the possibility of a metastable top squark. It is 
well known that the top squark can be the lightest supersymmetric particle 
(LSP) if the left-right stop mixing is comparable to the squark masses (for 



2 



reviews, see |3|). Combined with the requirement that the interactions be 
strictly R-parity invariant, this results in the stability of the top squark as 
the stop decays into other superparticles are not allowed kinematically. 

It is, however, also well known that such a scenario encounters serious 
cosmological problems. A stable stop would form stable hadrons. The abun- 
dance of such hadrons can be estimated in the standard cosmological model 
knowing the masses and cross sections of these particles. One finds that 
stable 100 GeV hadrons should exist in concentrations of one in 10 10 — 10 12 
nucleons ||. Thus they could be detected as anomalously heavy protons or, 
if they get bound in nuclei, they could be found in anomalous isotopes [|]. 
However, years of experimental searches for anomalous protons and isotopes 
have resulted in extremely tight bounds on their concentration in matter 



lj q=+l / Tlnucleon ^ 10 , 

<, q =-i/n nudeon < 4 x 10 -20 (1) 



for m=100 GeV; the displayed bound on n q= _i results from the search for 
anomalous isotopes of carbon. Apparently, these numbers are in conflict with 
the expected concentration of the stop containing hadrons. Thus, the top 
squark cannot be absolutely stable. 

The stability of the stop arises from the assumption of an exact R- 
invariance. However, recently, with the discovery of the neutrino mass, this 
assumption has lost some of its theoretical motivation. It seems more plausi- 
ble that R-parity is not an exact symmetry as is not lepton (baryon) number 
(see || for a review). In this case the LSP decays into ordinary particles. If 
R-parity is broken only slightly, the decay time can be very long. This may 
lead to other cosmological problems. For example, if the decay products 
contain quarks, the balance between protons and neutrons would be broken 
during the nucleosynthesis leading to unacceptable abundances of helium and 
other light elements |K| . This places strict constraints on the lifetime of the 



dark matter candidates such as neutralinos and gravitinos. These considera- 
tions can also lead to significant constraints on the lifetimes of new strongly 
interacting particles if one assumes that the baryon asymmetry in the "new" 
sector is the same as it is in the regular matter sector. This assumption, how- 
ever, is strongly model-dependent and needs not be the case. In particular, 
if the asymmetry in the stop sector is much smaller than that in the matter 
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sector, the concentration of the susy hadrons at the time of nucleosynthesis 
is too small (~ 10~ 10 ni, ar fj) to affect the formation of light elements. In 
this case the stop lifetime is left essentially unconstrained as long as 

T s top < 10 10 years . (2) 

The top squark with such a lifetime is the subject of our paper. 

The stop is allowed to decay into quarks and charged leptons via the R- 
breaking interactions: 

Wft = KjkUQA + KjADjDk , (3) 

where L i; Qi, Ui, and £>i are the usual MSSM superfields. The corresponding 
decay width is given by 

„ I A' , A'" 2 



-m; . (4) 
16tt * K J 

The requirement that all the stops have decayed by now can then be trans- 
lated into the lower bound on the couplings: 

A', A" > 1(T 21 (5) 

with mi = 100 GeV . The most recent upper bounds on these couplings from 



particle experiments can be found in [II 



A model containing such small numbers may seem ill-motivated. However, 
very small R-parity violation can naturally arise in some flavor models with a 
discrete symmetry jTJ] , JT3J . In particular, in the model of Ref . |T3|j , the (S3) 3 



symmetry prohibits renormalizable R-breaking operators. R parity arises as 
an accidental symmetry and can be broken by nonrenormalizable operators 
induced at the Plank scale. Then the effective low-energy trilinear operators 
get suppressed by (Mf/M P i) 2 and Mf/Mpi, where Mf is the flavor symmetry 
breaking scale. With a sufficiently low My, the R-breaking couplings in the 
range given by Eq.^Jcan be naturally produced. For example, the lower bound 
on the A', A" corresponds to Mf ~ 10 9 GeV if the R-breaking interactions 
are suppressed by Mp t . This provides a theoretical motivation for studying 
a possibility of a long-lived top squark. 
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Table 1: Lightest hadrons containing a top squark. 

3 Lightest Stop Containing Hadrons 

Since any observable object has to be a color singlet, a metastable top 
squark must be confined in hadrons. Within O(l0~ 23 )sec it will form a bound 
state with either an antiquark - mesonino, or a couple of quarks - sbaryon. 
Clearly, the lightest mesonino (sbaryon) is just as stable as the stop itself. 
Masses of the lightest shadrons can be approximated within 1 GeV by the 
stop mass which we will set to 100 GeV, however splittings between them 
will be very important for phenomenology. The natural candidates for the 
lightest hadron are listed in Table 1 (this is to be supplemented with the 
charge conjugate states). 

In the presence of isospin symmetry breaking, T° is not a mass eigenstate. 
Indeed, due to the diagram in Fig.l, T° mixes with T° forming Majorana 
fermions T\ and T-i- 

Let us estimate how large is the mixing. The leading contribution is given 
by the gluino diagram in Fig.l. Further, in the limit m~ 3> m|, the mixing 
is dominated by the diagram with a gluino mass insertion (Fig.2). It is also 
reasonable to assume that the lightest stop is predominantly right handed. 
Then, the generated effective Hamiltonian is given by 

n eff « 9 ^(u R i R y , (6) 

where V is the mixing matrix at the gluino vertex. To convert this interaction 
into a mass term, one needs to know the matrix element (T°\H e ff\T°). In the 
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case of the Standard Model K — K mixing, the needed matrix elements can 
be derived from the PCAC relation and the vacuum saturation assumption. 
In our case we will simply parametrize the matrix element assuming 

(T°\(u R t R ) 2 \f°) = F 2 , (7) 

where F has dimension of mass[]. The size of F can be estimated by the 
characteristic hadronic matrix element of a heavy quark meson 
(9(1) GeV since this combination is constant in the heavy quark limit. 
As a result, we get the following mass matrix in the (T°,T°) basis: 

where e = g m '" F . In principle, V tu may contain a complex phase, leading to 
CP violation in mixing. However, in this paper, we restrict ourselves to the 
CP conserving case. The consequent mass eigenstates are maximal mixtures 
ofT°andf°: 

Ti = -^(T°-f°), mi=m -e, 

T 2 = -±={T° + f°) , m 2 = m + e . (9) 
v2 

It is interesting to compare the mass splitting between the two eigenstates 
to a typical mass splitting in the neutral meson systems. For example, in 
the Standard Model, the B° - B° mixing arises only through a box diagram 
resulting in (Am/mjeo ~ 1CT 13 . In contrast, the T° — T° mixing exists 
already at the tree level and 

*2 = l^JL ~ !0-» (10) 

m 6 rng-mi 

for V tu ~ V£ KM , m- g ~ 300 GeV, m~ t ~ 100GW and F ~ 1 GeV. The 
result, of course, is sensitive to F, but it is clear that Am/m is enhanced by 
orders of magnitude over that in the Standard Model. 

Let us now turn to the discussion of the isospin partner of T° - T + . Since 
the top squark can be treated as a spectator, the mass splitting in this 

x In our convention (T°|T°) = <5 (3) (p -p'). 
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isospin doublet can be attributed to the — m u mass difference and the 
corresponding Coulomb interactions. As a result, T + is unstable and j3- 
decays: 

T+ -> T° + e+i/ (11) 

(It is understood that T + actually decays into a linear combination of the 
mass eigenstates T\ and T 2 .) Clearly, this process is governed by the usual 
Fermi interaction in analogy with a neutron /3-decay: 

H = J= T + 7^(l - a 75 )T° e 7M (l - 75)^ , (12) 

where a ~ 1 accounts for hadronic effects. The decay width is given by that 
of a neutron /?-decay: 

r ~ ^ (13) 

60vr 3 V ; 

with A = m T + —m T o ~ m D + — m D o ~ 5MeV. The corresponding decay time 
is about 3 sec. So within a few seconds all of the T + 's will decay into the 
T°'sg. 

One can go further in eliminating candidates for the lightest shadron. Let 
us consider the mass splitting between the vector isotriplet V and the scalar 
isosinglet S + . The mass difference arises due to the spin flip of one of the 
quarks. This effect is accounted for in the quark potential model via the 



spin-spin interaction (see, for example, [15|). Such an interaction is part of 
the one gluon exchange potential and is responsible for the p — ir and A — N 
mass differences. Again, treating a spin stop as a spectator, we can write 
the hyperfine interaction as follows: 

H h . f . = k qq sl-s 2 5 {z \r) . (14) 

The constant k depends on the quark masses and on whether two quarks or 
a quark and antiquark are interacting, k qq = \k qq . The scalar product s~! • s~2 
is 1/4 for the vector state and —3/4 for the scalar state. This situation is 
very similar to what we encounter in the p — 7r system up to the substitution 
kqq — > k qq . Therefore, 

my — m s + ~ 2^ mp ~ m " K ^ ~ MeV . (15) 



2 This estimate is sensitive to the exact value of the mass difference. The actual decay 
time may be an order of magnitude larger. See |l4| for similar considerations. 
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Another way to estimate the V — S + mass difference is to exploit the heavy 



quark symmetry of the strong interactions fll6| . In the heavy quark limit, the 
interactions of a heavy quark are independent of its spin and flavor. Thus 
the heavy quark spectroscopy depends on the light quark configuration only. 
Using this fact, one can approximate the V — S + mass difference by that of 
the E and A baryons containing a heavy quark. This gives 

m v - m s + ~ 200 MeV , (16) 

which is about 30 % away from the previous estimate^. 

In both cases the mass splitting is well above the pion threshold, so V will 
decay into S + within about 10~ 23 sec: 

V -> S + + 7i (17) 

Thus we are left with only two potentially long-lived shadrons - T° and S + . 
In principle, there is a possibility that one of them will decay into another. 
Since they carry different baryon numbers, a proton must be emitted in such 
a decay: S + — > T° + p. It is, however, quite unlikely that the S + — T° 
mass difference is above 1 GeV. To see this, we can again exploit the spin 
independence of the heavy quark interactions and use an analogy with the 
heavy quark hadrons. The S + — T° mass splitting can be compared to that 
for the A6( c ) baryons and B(D) mesons. This estimate yields 

m s + - m T o ~ 400 MeV . (18) 

Thus, the decay S + —>■ T° + p is not allowed kinematically. 

To summarize, we have argued that there are four long-lived shadrons 
(including the charge conjugate states): and S . Their lifetimes are 
given by the stop lifetime. Other light shadrons decay into these states 
either strongly or via /5-decay. 



4 Possibility of Detection 

Supersymmetric particles have been the subject of collider searches for 
many years (for recent reviews see [T7\ and references therein). In most of 

3 To obtain this number, we use the mass of A c and the average mass of E c , spin 1/2 
and E c , spin 3/2 (which should be the same in the heavy quark limit) plfl , Anologous 
information for the 6-baryons currently is not available. 
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these searches it was assumed that the LSP is a neutral particle which eludes 
detection. This assumption leads to the celebrated signature of supersym- 
metry - missing transverse momentum and energy. The recent LEP and 
Tevatron bounds on masses of the supersymmetric particles can be found in 

US- 



In the scenario which we are discussing, the LSP is a charged strongly 
interacting particle. If it were sufficiently long-lived, it would leave a track in 
the tracking chambers and thus could be detected directly. In particular, if 
a heavy charged particle is produced at the Tevatron, it can be identified in 
the CDF detector due to its relatively low velocity and, consequently, large 
ionization deposition, dE/dx. CDF measures dE/dx in the silicon vertex 
detector as well as in the central tracking chamber. In Run I, no events in 
excess of the expected background have been found, setting a limit of 85 GeV 
on the mass of a unit-charged scalar triplet [Hi]. A more recent analysis PU 



increases this bound up to about 100 GeV 0. 

In the previous subsection we have argued that there are four metastable 
shadrons. This is the case when these shadrons are in isolation. However, 
when interactions with matter are turned on, only one of them survives. Let 
us consider the states tu and tu propagating in matter. Apparently, tu has 
a larger interaction cross section since the annihilation channel is open: 

tu + p ^tud + vr° , 2vr , ... (19) 
tu + p -> tud + 7 , 27 , ... (20) 

As a result, the interaction eigenstates are different from the mass or flavor 
eigenstates. This is a well known phenomenon in K-meson physics: the same 
effect is used to regenerate the K$ component via interactions with matter. 
The effective mass matrix in the (T°,T°) basis now takes the form 

M'=( m £ V P"/ 7 °\ = ( m -^ , (21) 

\e m J \ 1/ \e — ^7 m ) 

where we have taken into account only the largest absorptive effect. The 
constant 7 > depends on the material and energy of the shadron. The size 

4 Many of these bounds are based on signatures specific to certain models and do not 
hold in general. For example, the Tevatron bound on the stop mass is based on the 
decay t — > cx°, which cannot occur in our model. 

interactions with the material of the detector may weaken this bound by 30-40 GeV 

h. 
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of 7 can be estimated very roughly using the nuclear interaction length Xj. 
For protons propagating in iron, Xj ~ 17 cm [^TJ. This number triples for the 
T° as it contains only one light (anti)quark ||22j] , yielding 7717 ~ 10~ 15 GeV . 
Therefore, m^ <C e and the interaction eigenstates can be well approximated 
perturbatively in 7. The corresponding effective masses are 



m'i — (m — e) 



m' 9 = (m + e) 



1 _ !2 
2 

1 _ !2 
2 



(22) 



Apparently, both mesonino states eventually vanish in a medium converting 
into sbaryons. For example, in iron, the attenuation length would be about 
1 meter. At very low energies, the mesoninos may interact with matter 
primarily through the electromagnetic annihilation of u and u if the strong 
channel ([H]) is kinematically prohibited. 

Similar considerations can be applied to the sbaryon S~ which contains 
light antiquarks. As a result, the only stop-containing particle which does not 
get destroyed in matter is S + . All other shadrons propagating in a medium 
will go through a series of reactions and eventually end up as S + . 

Nevertheless, if the A' and A" are anywhere within 10 orders of magnitude 
from the lower bound (||), all of the charged shadrons but V (which decays 
too quickly) would leave tracks in the tracking chambers. Since they are 
expected to be very massive, they would pass through the electromagnetic 
calorimeter undetected. With some energy loss (which may be considerable 
for slow particles p2]| ) in the hadron calorimeter, a shadron will penetrate 
the detector and will be triggered on as a muon. 

Let us consider the interaction of shadrons with the hadron calorimeter 
in more detail. In collisions with the nucleons, the momentum transfer at 
typical Tevatron energies is very small and the energy available in the center 
of mass system is likely to be below the pion threshold [[TJ|. Thus, the 
inelastic scattering would primarily go through the annihilation of the light 
quarks and antiquarks as in Eq.|T9|. Since the S + — T° mass difference is 
estimated to be around 400 MeV, a few soft pions (or an equivalent amount 
of radiation) can be produced in the collision of T 12 and p. Note that, at low 
energies, only the shadron states involving light antiquarks (T + , S~ , etc.) 
would deposit any tangible amount of energy in the hadron calorimeter. 
However, at the Tevatron, the deposited energy is unlikely to exceed a few 
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GeV in any case and thus the signal would be hard to distinguish from that 
generated by a charged colorless LSP such as a stau, especially taking into 
account that the energy resolution is only 80 — 100%/v^E |fTT| . 



The situation will improve at the LHC The energy deposition in the 
calorimeter will be sufficient to distinguish between a strongly and a weakly 
interacting LSP. The energy loss in a material as a function of the mass and 



energy of a strongly interacting particle was studied in Ref.p3| based on 
kinematical considerations. For instance, it was shown that a 100 GeV par- 
ticle at E — 1 TeV would trigger a few 10 GeV hadron showers. Also, it was 
noted that, due to a slower velocity, the opening angle of the shower would 
be significantly wider than that of the shower induced by regular hadrons. 
These calorimeter data combined with the tracking (and timing) information 
would allow to eliminate the SM background and determine the nature of a 
new (meta)stable particle. 

Let us now briefly discuss other implications of the scenario. The stop may 
be sufficiently long-lived as to be treated as stable. One may ask whether 
there could be non-collider experiments which can detect shadrons produced 
in the collisions. In principle, this can be done via mass spectroscopy. If a 
produced shadron is sufficiently slow, it loses most of its energy through 



ionization and thus can be trapped in the metal plates of the detector [22 
eventually ending up as an S + . The isosinglet S + is very unlikely to be bound 
in nuclei due to the Coulomb repulsion and absence of the one-pion-exchange 
potential. Therefore, it should be found in isolation. Having picked up an 
electron, it will form an anomalously heavy hydrogen atom. If such atoms 
existed in an appreciable concentration, they would be detectable by mass 
spectrometry [[/J as they possess an anomalously small charge/mass ratio. 
We would like to conclude with a remark: 

"Curiosity is a delicate little plant which, aside from 
stimulation, stands mainly in need of freedom." 

— Albert Einstein 

We have benefited from very helpful correspondence with T. Cohen, H. 
Dreiner, X. Ji, A. Kronfeld, D. Stuart as well as from discussions with mem- 
bers of Virginia Tech and U. of Maryland high energy groups. We are par- 
ticularly thankful to Lay N. Chang, S. Eno, R.N.Mohapatra and G. Snow 
for enlightening discussions. 
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Figure 2: The dominant contribution to T° — T° mixing. 
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